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We investigate arrays of three traps with two fermionic or bosonic atoms. The tunneling inter- 
action between neighboring sites is used to prepare multi-site dark states for the empty site, i.e., 
the hole, allowing for the coherent manipulation of its external degrees of freedom. By means of 
an ah initio integration of the Schrodinger equation, we investigate the adiabatic transport of a 
hole between the two extreme traps of a triple-well potential. Furthermore, a quantum-trajectory 
approach based on the de Broglie-Bohm formulation of quantum mechanics is used to get physical 
insight into the transport process. Finally, we discuss the use of the hole for the construction of a 
coherent single hole diode and a coherent single hole transistor. 



I. INTRODUCTION 

The beginning of the 21st century has brought the de- 
velopment of techniques to isolate and manipulate indi- 
vidual neutral atoms [l] , allowing to follow a bottom-top 
approach where quantum systems acquire classical fea- 
tures as their size and/or their coupling with the envi- 
ronment increases. In fact, in the last few years a lot 
of attention has been devoted to the field of atomtronics 
[2] , where atomic matter waves in optical [SHS] , magnetic 
[6^ and electric [T potentials play an analogous role to 
electrons in electronic devices. Atomtronics has the im- 
portant advantage over electronics that neutral atoms are 
comparatively less sensitive to decoherence than charged 
particles, i.e., interaction with the "classical" environ- 
ment can be almost completely inhibited for the former. 
In fact, neutral atom devices based on the coherent tun- 
neling of matter waves are, in general, designed to take 
profit of their inherent quantum features. In this context, 
several proposals on diode- and transistor-like behaviors 
have been deeply investigated for ultracold atoms [8] and 
Bose-Einstein condensates [9^ both in double and triple 
well optical potentials as well as in optical lattices with 
applications ranging from atomic cooling to quantum in- 
formation processing. 

Although atomtronic devices have yet to be realized 
experimentally, the time for coherent atomtronics is al- 
ready here mainly due to the fact that techniques for 
cooling and trapping atoms are by now very well estab- 
lished [TJ [SHS]. Neutral atoms can be stored and ma- 
nipulated in optical lattices, standard dipole traps, and 
microtraps. In particular, magnetic and optical micro- 
traps offer an interesting perspective for storing and ma- 
nipulating arrays of atoms with the eventual possibility 
to scale, parallelize, and miniaturize the atomtronic de- 
vices. Moreover, optical microtraps can take advantage 
of the fact that most of the current techniques used in 
atom optics and laser cooling are based on the optical 
manipulation of atoms. In fact, the possibility to store 



and to selectively address single optical microtraps, as 
well as initializing and reading out the quantum states 
in each of the sites has been already reported [3^ 4 . 

In this context, there is a need for the development of 
novel techniques to control the coherent flow of matter 
waves in optical and magnetic traps based on tunneling 
devices. Recently, we introduced a set of coherent tools 
[10 , named three- level atom optics (TLAO) techniques, 
to efficiently transport matter waves between the two ex- 
treme traps of a triple-well potential via the tunneling 
interaction. This adiabatic transport process is the mat- 
ter wave analog of the very well known quantum optical 
STIRAP technique pjj and it is based on adiabatically 
following an energy eigenstate of the system, the so-called 
spatial dark state, that, ideally, only involves the vibra- 
tional ground states of the two extreme wells. Extensions 
of these TLAO techniques to atomic wave packets prop- 
agating in dipole waveguides [12], to Bose-Einstein con- 
densates [13 , to the transport of electrons in quantum 
dot systems [14 , and to superconductors [15 , have been 
later performed. Even very recently, by exploiting the 
wave analogies between classical and quantum systems, 
de Longhi et al. [16 have experimentally reported, for 
the first time, light transfer in an engineered triple- well 
optical waveguide by means of the classical analog of the 
matter wave STIRAP. 

In the first part of this paper we will consider an ar- 
ray of three traps with two neutral atoms and extend 
the TLAO techniques to the transport of the empty site: 
the hole, see Fig. [ija). In the second part of the paper, 
we will discuss the use of the hole as an active player in 
coherent atomtronic devices [2 . We will design a single 
hole diode, see Fig. TJb), by tuning the interaction be- 
tween the atoms, allowing for the hole transport to be 
successful in one direction but inhibited in the opposite. 
Furthermore, we will engineer a single hole transistor, see 
Fig. [TJc), by addressing and manipulating the spin of an 
individual atom and changing the symmetry of the two 
atoms spin state. 



The article is organized as follows. In Section [TT] we 
will introduce the physical system under investigation: 
an array of three traps with two neutral atoms. In Sec- 
tion III we will describe the adiabatic transport of a hole 



between the two extreme traps introducing the concept 
of a spatial dark state for the hole. Through a numerical 
integration of the Schrodinger equation, we will address 
the adiabatic dynamics of a single hole in a triple well 
potential. By means of a quantum-trajectory approach 
based on the de Broglie-Bohm formulation of quantum 
mechanics, we will elucidate some specific features of the 
adiabatic transport technique here proposed. This Sec- 
tion will be followed by a discussion. Section \P7\ on the 
use of the hole as a key element in the building up of 
a coherent single hole diode and a coherent single hole 
transistor. In Section |V| the concluding remarks are pre- 
sented. 



II. PHYSICAL SYSTEM 

The system under investigation is sketched in Fig.jlja) 
and consists of an array of three optical traps that has 
been loaded with two fermionic or bosonic atoms in the 
lowest vibrational levels. For the physical implementa- 
tion of the dipole traps, we will consider the 2D array of 
optical microtraps discussed in Refs. [3, 4 where single 
site addressing as well as the ability to approach columns 
(or, alternatively, rows) of traps yielding coherent atomic 
transport has been demonstrated [4 . Note that the use of 
the 2D array allows for the performance of the techniques 
here described in a parallel fashion, realizing multiple ex- 
periments simultaneously. 

We assume that tunneling between sites occurs only 
in the column movement direction (namely, x) and thus, 
the main dynamics will be restricted to one dimension. 
For two identical atoms of mass m, the dynamics of the 
system are governed by the Hamiltonian: 
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with V being the trapping potential and U the interac- 
tion between the two atoms. 

For simplicity, we assume that the three-trap array po- 
tential consists of truncated harmonic wells centered at 
positions {xq.} with i = 1, 2, 3: 

V{x) = ^mul mm{{x - xqJ^} , (2) 

where uj^ is the longitudinal trapping frequency of each 
trap. For the one-dimensional model to be valid, we as- 
sume a tight transverse confinement such that transverse 
excitations can be neglected, i.e., ujp ^ uOx^ being uop the 
transverse trapping frequency. In this one-dimensional 
model, the cold collisional interaction between the atoms 
can be modeled by a contact potential of the form [17 : 

U{xi,X2) = 2hasUJpS{xi - X2), (3) 

with as the s-wave scattering length. 
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FIG. 1. (color online) (a) Simplified three- level description 
of the system under investigation: a three-trap array loaded 
with two atoms, e.g., two fermions, and one empty site, i.e., 
a hole. Matter wave STIRAP techniques will be applied to 
transfer the hole between the traps. ) is the localized state 
for the (fermionic) hole being at trap i = 1, 2, 3 and Ji is the 
hole tunneling rate between traps i and i + 1. (b) Single 
hole diode: for identical trap-approaching schemes but with 
appropriately tuned interactions between two bosonic atoms, 
the hole will be transported from the left to the right trap, 
but the inverse process will not succeed, (c) Single hole tran- 
sistor for two non-interacting fermionic atoms: by performing 
identical trap-approaching schemes, the hole transport from 
left to right will succeed or will be inhibited depending on 
whether the spin state of the two atoms is symmetric or an- 
tisymmetric. 



III. ADIABATIC TRANSPORT OF HOLES 

Following standard ideas of solid state physics we ad- 
dress the present problem in terms of a hole for either 
fermionic or bosonic atoms [18 . For the hole description 
to be valid, the following conditions must be fulfilled: (i) 
each trap contains, at most, one atom; (ii) all atoms are 
cooled down to the vibrational ground state of each trap; 
and (iii) tunneling is adiabatically controlled to strongly 
suppress the probability of double occupancy. To satisfy 
condition (iii) for identical fermions, we assume that the 
two atoms have parallel spins such that the Pauli exclu- 
sion principle applies and double occupancy in the same 
vibrational state is strictly forbidden. For bosons, we re- 
quire for the 5-wave scattering length to be large enough 
to, in the adiabatic limit, inhibit double occupancy. 

For the two fermions case, the entire two-atom state 
must be antisymmetric which means that if their spin 
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state is symmetric (antisymmetric) then the spatial wave- 
function must be antisymmetric (symmetric). For the 
two bosons case, since the entire two-atom state must be 
symmetric, then their spin state and their spatial wave- 
function must have the same symmetry. Taking into ac- 
count that the dynamics we will simulate only involve 
the spatial wavefunction, we will distinguish between the 
two cases where: (i) the spatial wavefunction is antisym- 
metric (fermions with symmetric spin state or bosons 
with antisymmetric spin state), and refer to it as the 
fermionic case; and (ii) the spatial wavefunction is sym- 
metric (bosons with symmetric spin state or fermions 
with antisymmetric spin state), and refer to it as the 
bosonic case. In the fermionic case, the contact potential 
from Eq. ([s]) will not play any role in the dynamics since 
at xi = X2 the (antisymmetric) spatial wavefunction van- 
ishes. Furthermore, the fermionic and hardcore (strongly 
interacting) bosonic cases will present equivalent dynam- 
ics [19 due to the fermionic exchange interaction. 
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FIG. 2. (color online) Coherent transport of a fermionic hole 
in a triple well potential: (a) Temporal variation of the trap 
distances in terms of a (its inverse is the width of the ground 
state of each isolated trap) and ujx (the longitudinal trapping 
frequency), (b) Snapshots of the two-fermion joint probability 
distribution for the particular times indicated by the arrows 
in (a). Initial state: (/)f (xi,X2) = (xi,X2|^f). Parameter 
values: admax — 9, admin — 1.5 and uJxtdeiay — 120. 



Our first goal will consist in developing an efficient and 
robust method to adiabatically transport holes between 
the two extreme traps of the triple well potential, by 
applying the matter wave analogue [10] of the quantum 
optical STIRAP technique [11]. Let us consider, thus. 



three in-line traps with one empty site and two identi- 
cal fermions, each one in the vibrational ground state 
of the two remaining traps as shown in Fig. [ija) (the 
bosonic case will be discussed later on). In a three-state 
approximation, where the Hilbert space of the system is 
restricted to the lowest three energy eigenstates, the spa- 
tial wavefunction of the two atoms can be expressed in 
the following basis: 
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Ifi) = [|</'l)l|</'2)2 - |</'2>l|</'l>2] , 



(4) 
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where \(pj)k = \^j{x^t))k is the time-dependent state 
of the k-th atom localized in the j-th trap. States 
— l^f ^2, ^)) with i = 1,2,3 account for the 
fermionic hole being in the left, middle, and right trap, re- 
spectively. Note that Eqs. (|4j|6| are antisymmetric since 
we have assumed that the spin state of the atoms is sym- 
metric. For truncated harmonic traps, see Eq. ([2|, the 
Ji tunneling rate of a single atom between the ground 
states of two adjacent traps i and i + 1 is given by [TO] : 



Ji{adi) = uj^ 



-1 + e("^^)' [1 + adi (1 - eifiadi))] 
(e2M.)2 _ 1) /^2adi) 



where di = |xo.^i — . | and a = ^/muJx/fl. Therefore, 
in the hole basis {|^f ), l^f^), )} the dynamics of the 
system are governed by the Hamiltonian: 



^0 Ji{t) y 

TRAPS = ^ I Jl{t) y J2{t) 
J2{t) 



(8) 



being Ji{= Ji) [20^ the hole tunneling rate between two 
adjacent traps, see Fig.[TJa). One of the three eigenstates 
of Hamiltonian ([8| is the so-called spatial dark state [10] 
that only involves the two states where the hole is in the 
extreme traps: 

\D^m))) = cose(t) l^f ) - sine(i) 1^3^), (9) 

with tan0(t) = Ji{t)/ J2{t). The transport of the hole 
between the two extreme traps of the triple well potential 
consists in adiabatically following state \D) from ) 
to 1^3^) by smoothly varying the mixing angle from 
6 = to 6 = 7r/2. As in standard optical STIRAP 
[TT] . it is convenient to establish a general adiabaticity 
criterion given by Jmax^deiay > A with J^^^ = (Ji)max + 
(-^2) max ^ being a dimensionless constant that for 
optimized temporal delays and tunneling profiles takes 
values around 10. The generalization of Hamiltonian (|8| 
and the adiabatic transport process to a trap array of 
arbitrary length can be found in Appendix A. 

The previously discussed three-level approach has been 
introduced to illustrate the main ideas behind the hole 



transport. Nevertheless, and in order to be accurate, in 
what fohows we wih numericahy solve the Schrodinger 
equation in real space. Figs. [2ja) and |2|b) show an 
exact simulation, i.e., an ab initio numerical integra- 
tion of the Schrodinger equation with the Hamiltonian 
given in Eq. ([T]), of the adiabatic transport process of a 
single hole in a triple well potential with two identical 
fermions. The initial state is ^f(xi,X2) = (xi,X2|(^f) 
with adi = ad2 = 9 while the expected final state 
is, up to a global phase, (xi,X2) = (xi,X2|^f) with 
adi = ad2 = 9. For the truncated harmonic potentials 
here considered, at adi = 9 the tunneling rate between 
adjacent traps is almost negligible and they can be con- 
sidered as isolated. For the time variation of the trapping 
potential, we have taken the middle trap to be static at 
X = while displacing only the two extreme traps. Note 
that the hole transport sequence, Fig.[2|a), starts by first 
approaching the two occupied traps and later approach- 
ing the empty trap to the middle one. Fig. [2|b) shows 
different snapshots for the temporal evolution of the two- 
fermion joint probability distribution | ^(xi , X2 , | ^ . Note 
that the diagonal Xi = X2 is forbidden due to the Pauli 
principle and the probability density is mirrored at both 
sides of this diagonal due to the antisymmetrization of 
the wavefunction. 

As indicated in Eq. (|9| , the hole is transferred from the 
left to the right trap with an ideally negligible popula- 
tion in the middle one. Therefore, the signature that the 
hole has been transferred through the matter wave STI- 
RAP technique is that the counter-diagonal xi = —X2 is 
practically not populated (see the 3rd, 4th, and 5th snap- 
shots of Fig. [2|b)). However, resorting to the continuity 
equation associated with the two-atom matter wave, the 
corresponding wavefunction must cross at some point the 
forbidden counter-diagonal. To get physical insight into 
this particular feature of the adiabatic transport process, 
we will discuss now very briefiy the previous simulations 
by means of quantum trajectories d la de Broglie-Bohm 
[2T] . See Appendix B for details of the quantum trajec- 
tories formulation. 

Fig. [3] shows a set of quantum trajectories calculated 
from the time evolution of Fig. [2|b). Their initial posi- 
tions, see Fig. p^a), are randomly distributed according 

to l^f (xi,X2)| . As expected. Fig. jsj^b) reveals that the 
time evolution of the quantum trajectories follows the 
evolution of the wavefunction, ending up distributed ac- 
cording to |(/)3^(xi, X2) 1^, see Fig. [sj^c). In Fig. [sj^b) we 
have also plotted the atomic probability distribution for 
the intermediate time Uxt = 360 corresponding to the 4th 
snapshot in Fig. [2|b). Clearly, when crossing the forbid- 
den counter-diagonal, each quantum trajectory suddenly 
increases its velocity, see Fig. [sj^d), such that the den- 
sity of trajectories per unit time vanishes in this counter- 
diagonal. In addition, quantum trajectories make a 
detour from the central region of the counter-diagonal 
where the probability distribution is significantly smaller. 

For the characterization of the transport probabilities 
between the different traps, after applying the matter 
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FIG. 3. (color online) de Broglie-Bohm trajectories corre- 
sponding to the temporal evolution of the system shown in 
Fig.l2|b). (a) Initial distribution of quantum trajectories (cf. 
Fig. |2[b-l)). (b) Evolution of the trajectories in configura- 
tion space together with the joint probability distribution at 
(jJxt — 360. (c) Final distribution of trajectories (cf. Fig.[2|b- 
6)). (d) Evolution of the trajectories as a function of time. 
In order to allow for the easy visualization of the transport 
process, (b) and (d) only show a reduced number of quantum 
trajectories. 



wave STIRAP sequence shown in Fig. |2ja), we define 
here their associated fidelities. These definitions will be 
also very useful for the characterization of the atomtronic 
devices introduced in Sec. HVl 

Thus, by starting with the hole in trap i = 1, 2, 3 and 
performing the temporal evolution, we check the popu- 
lation of each trap j = 1,2,3. We denote, for either a 
fermionic (F) or a bosonic {B) hole case, the state of the 
system at the end of the process as ^^(xi, X2, T), be- 
ing T the total duration of the STIRAP sequence. The 
population in the localized state (j)^^^ {xi^X2) after this 
evolution will be given by the product between this state 
and the evolved one, namely: 



F/B 



1^3 



/>f^^ {xi,X2)(t)f^'^ {xi,X2,T)dxidx2 



Thus, R 



F/B 



(10) 

is the fidelity of the fermionic (F) or bosonic 
(B) transport process of the hole between traps i and j 
after applying the hole matter wave STIRAP sequence, 
see Fig. |2|a). Note that for the matter wave adiabatic 
transport process, we want to maximize F^_^. 

Fig. |4|a) depicts F[_^^ in the parameter plane 
{t delay, dmin} for the holc transport process (see Fig. ma) 
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FIG. 4. (color online) Matter wave STIRAP fidelity: (a) Ff_^s 
in the parameter plane {t delay, dmin}- (b) Ff_^^ in the param- 
eter plane {as, dmin} for as positive and (c) for as negative. 
Scattering lengths of ^^Rb {aas = -7.98 • 10~^) and ^'^Rb 
{aas = 2.32 • 10~^) are indicated by dashed black lines. The 
temporal variation of the traps and the rest of the param- 
eters is given in Fig. |2ja). Transverse trapping frequency 
(jjp = 24(jJx 4 ^^^5 bosonic case, uJxtdeiay = 120. The 

bright (yellow) area indicates the area where the fidelity is 
larger than 0.99. 



for the definition of these two parameters). It comes clear 
that for a large set of parameters the fidelity is larger 
than 0.99 (see the bright yellow area in Fig. [4]^ a)), show- 
ing that the hole transport from left to right via matter 
wave STIRAP is a robust and efficient technique, pro- 
vided that the adiabaticity condition is fulfilled. 

We have also simulated the hole transfer process for 
the case of two bosonic atoms by integrating the corre- 



sponding Schrodinger equation. In this case, the local- 
ized states for the bosonic hole, l^f ), are given by the 
symmetrized versions of Eqs. (|4]-[6|. Figs. |4jb-c) show 
Fi_^2, parameter plane {a dmin}- As shown in 

the figures, the adiabatic transfer process succeeds for 
aas = —7.98 • 10~^ and aas = 2.32 • 10~^ corresponding, 
respectively, to the s-wave scattering length of ^^Rb and 
^^Rb [3j, while it breaks down for weaker interactions 
since then double occupancy starts to play a dominant 
role. For large absolute values of the 5- wave scattering 
length, bosons become hardcore and then their dynamics 
is equivalent to that of the fermionic case. 



IV. ATOMTRONICS WITH HOLES 

Making use of the fact that the hole transfer process 
here presented is spatially non- symmetric, we will now 
discuss both a coherent single hole diode and a coherent 
single hole transistor in a triple-well potential. The hole 
transfer from left to right and vice versa strongly depends 
on both the two-atom collisional interaction and the ex- 
change interaction and, therefore, both interactions will 
be used here to control the diode and transistor operation 
regimes. 



A. Single hole diode 

In this section we will design a single hole diode by 
using the collisional interaction between two bosons as a 
control parameter to allow the hole transport from left 
to right, and inhibit the transport from right to left, see 
Fig.jljb). Thus, Fig.jsja) shows the fidelity of the Dosonic 
hole transport processes Ff_^^, against 
the strength of the s-wave scattering length. The pa- 
rameter values for the temporal variation of the traps 
are taken as in Fig. [2] such that the fidelity of the hole 
transport process from left to right, ^1^3 (red circles in 
Fig. [sja)), is larger than 0.99 above a certain threshold 
value for the scattering length, indicated by point A in 
Fig. [sja), i.e., when the interaction is strong enough and 
the bosons become hardcore. 

By performing the same trap-approaching scheme but 
with the hole starting on the right trap, the process cor- 
responding to the hole being transferred from the right 
to the left trap (green triangles in Fig. p[a)) is inhibited 
{F^_^^ ~ at point B) or succeeds (i^^i 1 at point 
C), depending on the value of aas. 

We thus define the fidelity of the diode as: 

Fo = F,^_,{l-Fi_,). (11) 

since this fidelity is maximal when the bosonic hole is 
transported from left to right and, simultaneously, the 
opposite process consisting in the hole transport from 
right to left is inhibited, i.e., F]j = 1 when ^1^3 = 1 and 
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FIG. 5. (color online) Single hole diode: For a bosonic hole 
and as a function of the s-wave scattering length: (a) Fi- 
delities of hole transport processes Ff_^s (I'ed circles), Ff_^i 
(green triangles) a nd F ^_^2 (blue squares) and (b) diode fi- 



delity Fd (see Eq. (11) in the text). See text for the defini- 
tions of points A, B and C. The rest of the parameter values 



are as m 



As can be seen in Fig.jsjb), that shows Fjj as a function 
of the scattering length, by tuning the product of the in- 
verse of the size of the wavefunction and the 5- wave scat- 
tering length to aag ~ 4.25 • 10~^, where Fjj ~ 1 (corre- 
sponding to point B), we obtain a scheme that transports 
the hole from the left to the right trap but transfers a hole 
from the right trap to the middle one (see blue squares in 
Fig. [sja)). Note that an ideal diodic behavior where the 
hole ends at the right trap, no matter if initially it was 
at the left or the right trap, would violate the unitarity 
of the quantum evolution. 



B. Single hole transistor 

Fig. |4ja) shows that in the fermionic case, the hole 
transport from left to right achieves high fidelities. On 
the other hand. Figs. [4]^b-c) show that in the weakly in- 
teracting bosonic case, i.e., for — > 0, the hole transport 
does not perfectly succeed. In fact, we have checked that 
for as = the fidelity i^i^s vanishes with the hole ending 
in a superposition between being in the left and middle 
traps. From the previous observations, it is possible to 
design a single hole transistor where the spin state of the 



atoms is used to control the hole current from the left to 
the right trap. 

For instance, it is straightforward to check that, for 
two fermions in the middle and right traps (hole in the 
left), the state with symmetric spin state l^*) = (|t)i|t 
)2 — |i)iJi)2)/A/2 and antisymmetric spatial state (|^f ), 
cf. Eq. Q), 

1. 



1^)1^0 = 2[lt^2)i|t ^3)2 



It ^3)111^2)2 

-i;(/:'2)ii;(/:'3)2 + i;(/^3)ii;(/^2)2], (12) 

and the state with antisymmetric spin state \A) = (|i)i|t 
)2 — |t)i|i)2)/>/2 and symmetric spatial state (|^f )), 

l^)l^f ) = l[\i ^2)l|t ^3)2 + 1; (^3)l|t ^2)2 

-|t^2)i|;(/:'3)2-|t(/:'3)i|;^2)2], (13) 

are coupled via a spin flip on the atom in the middle trap, 
i.e.. It V^2)fe ^ |i ^2)/c- A similar argument for bosonic 
atoms can be done between states |*S')|^f ) and |A)|^f ). 

This control over the system behavior between the 
bosonic and fermionic cases allows us to create a coher- 
ent hole transistor scheme where the matter wave STI- 
RAP sequence from the left to the right trap succeeds 
or is inhibited depending on the spin state of the atoms. 
The case for two non-interacting fermions is depicted in 
Fig. [ijc). Thus, the fi gure of merit corresponds to the 
maximization of the transistor fidelity defined as: 



Ff_3 (1 - Ff_3) 



(14) 



since it will be maximal when Ff_^^ ~ 1 and Ff_^^ ~ 0. 
As we have discussed, for parameters of Fig. [2|a), Ft > 
0.99 for non-interacting atoms. 

To further test the fidelity of the atomtronic transistor, 

F IB 

we have calculated Ft, computing ^1^3 , adding a jitter 
in the trap positions (to simulate experimental imperfec- 
tions), see Fig. l6Fa), such as: 



d\(i) = d^{t) + As cos Ust, 



(15) 



being d^{t) the distance between the i-th and z + l-th trap 
shown in Fig.|2|a). As and ujs represent, respectively, the 
amplitude and the frequency of the jitter. The results 
of Ft for different values of As and Us are plotted in 
Fig. [oj^b) and indicate that for small jitter amplitudes 
the transistor still works with high fidelity {Ft > 0.99) 
for a wide range of frequencies, except for those close 
to Ux (trapping frequency) that, as expected, excite the 
atoms to unwanted vibrational states. 



V. CONCLUSIONS 

Using the matter wave analogue \iOl of the quantum 
optical STIRAP technique [TT], we have proposed an ef- 
ficient and robust method to coherently transport empty 
sites, i.e., holes, in arrays of three dipole traps with two 
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FIG. 6. (color online) Single hole transistor: (a) Tempo- 
ral variation of the trap distances shown in Fig. [2[a) with 
an added jittering of amplitude a As = 0.3 and frequency 
Us — O.lcJcc- The inset shows a sketch of the three traps and 
its simulated jittering, (b) Fidelity Ft of the transistor (see 
Eq. (14)) as a function of the amplitude of position jitter for 
different jitter frequencies, as = 0, and the rest of the param- 



eter values are as in Fig. |2] 



neutral atoms. The coherent transport process consists 
in adiabatically following a spatially delocalized dark 
state by an appropriate temporal control of the tunneling 
rates. We have first introduced the transport process in a 
simplified three-state model to, later on, simulate it with 
exact numerical integrations of the Schrodinger equation 
for time dependent potentials. Some particular features 
of the adiabatic matter wave dynamics, such as the trans- 
port through forbidden regions, have been elucidated by 
means of de Broglie-Bohm quantum trajectories. Finally, 
making use of both the collisional interaction and the ex- 
change interaction, we have analyzed in detail hole trans- 
port schemes for the implementation of a coherent single 
hole diode and a coherent single hole transistor in a triple 
well potential with two neutral atoms. 

Additionally, we want to note that it is possible to engi- 



neer a quantum information processing scheme that uses 
the hole as the qubit carrier, the computational states 
defined by the presence of the hole in one out of the 
two extreme traps of a triple well potential. The imple- 
mentation of single qubit gates could be performed by 
taking advantage of the tunneling between traps, and a 
controlled phase gate between two qubits could be imple- 
mented by using the collisional interaction between the 
atoms, in a similar manner as described in Ref. [22 . In 
fact, the hole transport process discussed along the pa- 
per could be used to coherently prepare particular qubit 
states for the hole. 
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Appendix A: Single hole in an array of n (odd) traps 

In most quantum computation proposals with trapped 
neutral atoms, a defect-free quantum system where all 
sites of the lattice are occupied by exactly one atom is 
needed to start the information processing, e.g., empty 
sites must be removed from the physical area of com- 
putation. With this purpose, here we will extend the 
adiabatic transport method presented in the main text 
to a trap array of arbitrary length by means of the Hub- 
bard model, generalizing Hamiltonian ([8| to a system of 
n (odd) traps. This technique could be used to transport 
holes away from the area of interest in order to prepare 
defect-free trap domains to, eventually, perform quantum 
computations. 

Note that the exact numerical simulation, i.e., by 
means of the corresponding Schrodinger equation, of the 
adiabatic transport of a hole in an array with a large 
number of single-occupancy traps is extremely demand- 
ing from a computational point of view and out of the 
scope of the present paper. However, the Hubbard for- 
mulation for the hole that we will now introduce could 
be used to simulate the hole dynamics following the lines 
of previous works on multi- level STIRAP [23] , 

Let us consider a single-occupancy ID array of n (odd) 
traps loaded with n — 1 spin-polarized fermions (or hard- 
core bosons) that, therefore, presents one isolated defect 
consisting of an empty site in one of its extremes. Again, 
the goal is to adiabatically transfer this empty site, i.e., 
the hole, from one extreme of the array to the other. The 
dynamics of a spin polarized non interacting Fermi gas 
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loaded in a ID trap array is described by the fermionic 
Hubbard Hamiltonian [24^ : 



(Al) 



Operators cj and q are the fermionic creation and anni- 
hilation operators at site i satisfying fermionic anticom- 
mutation relations and Ji accounts for the tunneling be- 
tween neighboring sites (see Eq. We have dropped 
the on-site energy term [24] since we are considering a 
homogeneous system with fixed number of atoms. 

In this context, we consider a hole (an empty site) as 
a virtual particle whose vacuum state, corresponds 
to all sites occupied with one fermion, reading: 



\n) = clcl..ci\n), 



(A2) 



with \Q) the fermionic vacuum state. Since we do not al- 
low for transitions to excited vibrational states and con- 
sider n traps with n — 1 atoms, the dynamics of our sys- 
tem is constrained to remain in states with a single hole, 
{Cj\Q)}, with cj = Ci being the hole creation operator at 
site i. Then, in terms of these on-site hole operators, and 
hole tunneling rates Ji{= Ji), Hamiltonian (Al) reads: 



H 



C]Ci. 



) . (A3) 



For n = 3 we retrieve Hamiltonian ([8|. For n odd, it 
is straightforward to check that this Hamiltonian has an 
energy eigenstate: 

_ ~2~ / m-1 ^ 

\D) = y: i-ir^' n -^2^-2,-1 

m=l = l 

n ^2m+2J CL-ll^), (A4) 

j=o J 

that satisfies H\D) =0. Note that \D) not only involves 
the first and last traps, but also all Cl\Q) with i odd. In 
this case, the hole transfer following state \D) would be 
achieved by favoring first the tunneling rates Ji with even 
i and then the ones with odd i [23i |25j . Finally, notice 



also that hole transport based on the adiaba tic fo llowing 
of the multi-site spatial dark state given in (A4) can be 
straightforwardly applied to the case of hardcore bosons. 



Appendix B: de Broglie-Bohm formulation 

By casting the polar form of the wavefunction, cp = 
^^^s/n ^ into the Schrodinger equation (with the conven- 
tional meaning of the symbols). 



..dip 



2m dx\ 2m 
and separating real and imaginary parts, we obtain 



(Bl) 



dS_ 
'dt 



' dt 



= V- 



E- 



2m V dxi 

i=l,2 

^ 2mR~dxf' 

i=l,2 * 



m oxi V (^Xi 

i=l,2 * ^ * 



(B2) 



(B3) 



Eq. (B2) is the so-called quantum Hamilton-Jacobi 



equation because of its similarity with the (classical) 
Hamilton-Jacobi equation but with one additional term, 
the quantum potential, accounting for the quantum fea- 
tures of the system. This similarity suggests the defini- 
tion of the particle velocity as: 



V^[t] = 



1 dS 

m dxi 



(B4) 



{xi[t\,X2[t\,t) 



Thus, Eq. (B3) becomes a continuity equation ensuring 



that the trajectories distribution i s giv en by i?^(xi, X2, t) 
at all times. After solving Eq. (Bl) (or, alternatively, 
Eqs. ( |B2| ) and ( |B3[ )) and distributing the initial positions 
of trajectories following the probability density function 
i?^(xi, X2, to), we find the quantum trajectories (time 
evolution of the positions) as: 



Xi[t] = / Vi[t]dt. 

J to 



(B5) 



We have followed the previous approach to obtain the 
quantum trajectories displayed in Fig. [3j 
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